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ABSTRACT 

The initial mass function (IMF) for the stars is usually fitted by three straight 
lines, which means seven parameters. The presence of brown dwarfs (BD) in- 
creases to four the straight lines and to nine the parameters. Another common 
fitting function is the lognormal distribution, which is characterized by two pa- 
rameters. This paper is devoted to demonstrating the advantage of introduc- 
ing a left truncated beta probability density function, which is characterized 
by four parameters. The constant of normalization, the mean, the mode and 
the distribution function are calculated for the left truncated beta distribution. 
The normal-beta (NB) distribution which results from convolving independent 
normally distributed and beta distributed components is also derived. The chi- 
square test and the K-S test are performed on a first sample of stars and BDs 
which belongs to the massive young cluster NGC 6611 and on a second sample 
which represents the star's masses of the cluster NGC 2362. 

Subject headings: Stars: luminosity function, mass function; Stars: fundamental 
parameters; Methods: statistical 



1. Introduction 

The distribution in mass of the stars has been fitted with a power law starting with 
Salpeter (1955). He suggested ^(m) oc where ^(m) represents the probability of having 
a mass between m and m + dm\ He found a = 2.35 in the range IOMq > M > IMq; this 
value has changed little with time and a recent evaluation quotes 2.3, see Kroupa (2001). 
Subsequent research has started to analyze the initial mass function (IMF) with three power 



- 2 - 



laws, see Scalo (1986); Kroupa et al. (1993); Binney & Merrifield (1998) and four power 
laws, see Kroupa (2001). A first comment on this temporal evolution is that the name 
is not appropriate because the power function distribution ^(m) cx is defined only for 
positive values of b, see Evans et al. (2000). A second comment is that the exact name for 
a probability density function (PDF) ^(m) oc m"'-'^"'"^^ with c > is the Pareto distribution. 
A third comment is that this progressive increase in the number of segments has limited 
the development of new or modified PDFs. The approach to the IMF by a continuous 
distribution has been modeled by the lognormal distribution in order to fit both the range of 
the stars and the brown dwarfs (BDs) regime, see Chabrier (2003). Recall that usually the 
standard PDFs such as the lognormal, gamma, generalized gamma, and WeibuU are defined 
in the interval < a; < oo. The fact that the number of stars with mass m < O.OTMq is 
nearly zero suggests a left truncated PDF. Our analysis has therefore been focused on the 
beta distribution, which by definition has an upper bound. From the previous analysis the 
following questions can be raised. 

• Is it possible to find the constant of normafization for a left truncated beta PDF? 

• Is it possible to derive an analytical expression for the mean, the mode, and the dis- 
tribution function of a left truncated beta PDF? 

• Is a left truncated beta PDF an acceptable model for the IMF as well as a real sample 
of masses? 

In order to answer the previous questions, we first review some standard PDFs, sec Section 
2. We subsequently introduce the various beta PDFs, the convolution of a beta PDF with 
a normal PDF, and a left truncated beta, sec Section 3. In order to find out which PDF 
performs best, the two main criteria which report the goodness of fit are found in Section 4. 
A comparison between various continuous PDFs and the left truncated beta is carried out 
in Sections 5.2 and 5.3 for two samples of stars. 

2. Distributions commonly used 

This section reviews some standard PDFs, namely, the lognormal, gamma, generalized 
gamma, Pareto, truncated Pareto, and the recently developed Double Pareto-lognormal 
distribution. 
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2.1. Lognormal distribution 

Let X be a random variable taking values x in the interval [0, oo]; the lognormal PDF , 
following Evans et al. (2000) or formula (14.2)' in Johnson et al. (1994), is 

/.« = ^exp^ii«! , (1) 

or 

1 - {In X - uln^ 
Iln = 7= exp — , 2 

where m — exp //^^v and ^iln — logm. 



2.2. Gamma distribution 

Let X be a random variable taking values x in the interval [0, oo]; the gamma PDF is 

p{x-h,c)^^^^L__ , (3) 
where T{z) is the gamma function 

poo 

V{z) = / e-H'-^dt , (4) 



see formula (17.1) in Johnson et al. (1994). 



2.3. Generalized gamma distribution 

Let X be a random variable taking values x in the interval [0, oo]; the generalized gamma 
PDF, following Evans et al. (2000), is 

f{x; a, b, c) = c ^^^^^^ x"'^ exp {-hx") , (5) 
see formula (17.116) in Johnson et al. (1994). 



2.4. The Pareto and the truncated Pareto distributions 

Let X be a random variable taking values x in the interval [a, oo], a > 0. The Pareto 
PDF is defined by 

f{x; a, c) = ca"x-("+^) , (6) 
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with c > 0, see formula (20.3) in Johnson et al. (1994). The traditional Salpeter slope is 
therefore -(c+1). An upper truncated Pareto random variable is defined in the interval [a, b] 
and the corresponding PDF, following 

Goldstein et al. (2004); Aban et al. (2006); Zaninetti & Ferraro (2008); White et al. 
(2008), is 



Their means are 



and 



frix; a, 6, c) = ^ _ . (7) 

E{x;a,c)^-^ (8) 
c — 1 

ca (-1 + (f)-') 



2.5. The double Pcireto-lognormal distribution 

The double Pareto lognormal distribution has been recently derived, see formula (22) 
in Reed & Jorgensen (2004), and has been used to fit the actual sizes of cities , see Giesen 
et al. (2010) ; its PDF is 

fix; a, ^, /X, a) = 1/2 (e^^^C-^+^M-^ H^))erfc(l/2 i^<^' + 1^ 'H^))^ ^ 

a 

(J 

where a and (3 are the Pareto coefficients for the upper and the lower tail, respectively, 
and (7 are the lognormal body parameters, and erfc is the complementary error function. 
The parameters can be found minimizing the maximum distance, D, of the K-S test, see 
Sect. 4. 



3. Various Beta distributions 

This section reviews the beta PDF defined in [0, 1], the beta with scale PDF defined in 
[0, h] and the general beta defined in [a, b]. The left truncated beta PDF defined in [a, b] but 
with a finite value of probability at x = a is explored. The convolution of a beta distribution 
with a normal distribution is also discussed. 
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3.1. Beta distribution 

Let X be a random variable taking values x in the interval [0, 1]; the beta PDF is 

with a > and /5 > 0, see Evans et al. (2000) or formula (25.2) in Johnson et al. (1995). 
Here B is the beta function defined by 

B(a, h) ^ f - tf-^dt = ^^^^^ . (12) 

Jo 



Its mean is 



and its variance, 



r(a + b) 

E{x;a,P)^^ , (13) 

a\x;a,l3)^ — ^ , (14) 

see formula (25.15a) in Johnson et al. (1995). The mode is at 

a - 1 

m{x-a,^)^- — — . (15) 

a — 2 + p 

The method of matching moments gives the following parameter estimation 

f X il — x) \ , ^ 

^-^"( ,2 -Ij ' (16) 



and 

P^(l-x) ('-^^^ 

where x and are the mean and the variance of the sample. The distribution function (DF) 
is 

DF(x-a 6)^^^^^^^^^^1^^±^^^±^ (18) 

where 2-^1(0, b; c; z) is the regularized hypergeometric function Abramowitz & Stegun (1965); 
von Seggern (1992); Thompson (1997); Gradshteyn, I. S. and Ryzhik, I. M. and Jeffrey, A. 
and Zwillinger, D. (2007); Olver et al. (2010). 



3.2. Beta distribution with scale 

Let X be a random variable taking values x in the interval [0,6]. The beta with scale 
PDF is 

Its expected mean is 

o. b 

E{x;b,a,/3)^^- , (20) 



a{x;b,a,P)l = - tL— . (21) 



and its variance, 

^2 /3ab'^ 
~ {l + a + /3){a + /3) 

The mode is at 

6 (a - 1) 



m{x-b,a,/3)t,^ o ■ (22) 

p — z + q; 

The DF is 

DFfc a;; 6, a, /3 = — — . (23) 

p (a, p) q; 

The three parameters can be estimated by 

b — maximum of sample , (24) 



X i-bx + a;^ + s^) 
« = ^l, ' (25) 



{b - x) {-bx + x^ + s^) 



(26) 



3.3. General Beta distribution 

Let X be a random variable taking values x in the interval [a, b]. The general beta PDF 

is 

/..(.;a,M,/^)^ (^-")(---Q^--^7 , (27) 
see formula (25.1) in Johnson et al. (1995). Its expected mean is 

E{x]a,b,a,^)ab ^ °^^'^^f > (28) 
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and its variance, 

a(a;;a,6,a,/5)^, = --^-— ^-^^ . (29) 

The mode is at 

m{x;a,b,a,P)ab^ — ■ (30) 

—2 + a + p 

The four parameters can be estimated by 

a = minimum of sample b = maximum of sample , (31) 

{—X + a) ( —X b + b a + x^ + s^ — d x] 
« = ^ 7 ^ , (32) 

S2 (d - bj 

lb — x] ( —X b + b d + x'^ + s'^ — d x) 

7 ^ , (33) 

S2 (d - b) 

see the discussion in section 25.4 of Johnson et aL (1995). 



3.4. Left truncated beta distribution with scale 

Let X be a random variable taking values x in the interval [a, b] and having a finite 
value in a. The left truncated beta with scale PDF is 

frix; a, b, a, /3) = K x""^ (6 - xf~'^ , (34) 

where the constant is 

^ -aT{a + f3) 

bt'-^H a"r {a + 13)- fe/^-i+ap (1 + a) T (/3) ' ^ ^ 

and 

i7 = 2Fi(«,-/3 + l; 1 + ^) . (36) 

The constant of normalization can be obtained from the integral of the beta with scale PDF 
as represented by Equation (19). The integral 3.194.1 on p315 of Gradshteyn, I. S. and 
Ryzhik, I. M. and Jeffrey, A. and Zwillinger, D. (2007), 



r x"-^ {1 + /3 xY dx 
Jo 

e/^i-H2Fi(/x,-z/; 1 + /x; -^u) 



(37) 
1^ 
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is useful in the analytical derivation of the constant. This PDF, which is new and therefore 
cannot be found in Johnson et al. (1995), at x = a is not zero but takes the finite value 

fT{a;a,b,a,l3) = Ka^'-^b-af^^ . (38) 

Its expected mean is 

y-^b'+'^2Fi{-/3 + 1, 1 + a; 2 + a; 1) 



E(x; a, b, a, P)t — K 



-fe/3-iai+«2Fi(-^ + 1, 1 + a; 2 + a; f ) 



1 + a 

The mode is at 



1 + a 

) ■ (39) 



m{x-b,a,(])T = -^^^^r^ , (40) 
a — Z + p 

and in order to exist one must have rriT > a. The DF is 

DFt{x; a, b, a, (3) = 

T{a + (5) a;"2Fi(a, -(3 + 1-1 + a- t)Hd'' 

b 

-2Fiia, -P + l;l + a; ^)a"r {a + /3) Hi,-' , (41) 

where 

2Fi{a,-p + l;l + a;^)a'^r{a + p)-r{p)bT{a)a . (42) 
The survival function is 

ST{x;a,b,a, /3) = 1 — DFT{x;a,b,a, (3) . (43) 

The four parameters can be estimated by 

d — minimum of sample b — maximum of sample . (44) 

A first couple for a and (3 can be obtained from those of the beta distribution with scale as 
given by Equations (25) and (26). A subsequent numerical loop around the previous values 
gives the couple which minimize the x^- 
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3.5. Beta distribution + normal 

We consider the sum Z = X + Y where X is a standard normal random variable, 
N{x;a), and F is a general beta distribution, fab{y',a,b,a, (3), as represented by the PDF 
(27). The sum, NB{z; a, b, a, (3, a), is 

NB{z;a,b,a,P,a)= f N{z - y)fabiy; a,b, a, /3)dy . (45) 

J a 

A similar example, uniform + normal, can be found in (Brandt & Gowan 1998, Sec. 5.11.2). 

This integral has an analytical solution for a and /3 integers, for example, when a=l 
and /3—1, 



NB{z;a,b,l,l,a) = 

1/2 erf I 1/2^-1/2^ | (a - 6)"' - 1/2 erf | 1/2^-1/2^) (a - 6)"' , (46) 



where erf is the error function. 



4. Goodness of fit tests 

The occasional reader may question which is the best fit for the distributions analyzed 
here. In order to answer this question, we first introduce x^, which is computed according 
to the formula 

X^ = E^^^T^, (47) 

where n is the number of bins, Tj is the theoretical value, and Oi is the experimental value 
represented by the frequencies. The theoretical frequency distribution is given by: 

Ti = NAxip{x) , (48) 

where N is the number of elements of the sample, Axj is the magnitude of the size interval, 
and p{x) is the PDF under examination. The size of the bins, Axj, is equal for each bin in 
the the case of linear histograms, but different for each bin when logarithmic histograms are 
considered. 

A reduced merit function Xred evaluated by 



(49) 
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where NF — n — k is the number of degrees of freedom, n is the number of bins, and k is 
the number of parameters. The goodness of the fit can be expressed by the probabihty Q, 
see equation 15.2.12 in Press et al. (1992), which involves the degrees of freedom and the 
X^- According to Press et al. (1992), the fit "may be acceptable" ii Q > 0.001. The Akaike 
information criterion (AIC), see Akaike (1974), is defined by 

AIC ^2k- 2ln{L) , (50) 

where L is the likelihood function and A- the number of free parameters in the model. We 
assume a Gaussian distribution for the errors and the likelihood function can be derived from 
the statistic L oc exp(— ^) where has been computed by Equation (47), see Liddle 
(2004), Godlowski & Szydowski (2005). Now the AIC becomes 

AIC ^2k + x^ ■ (51) 

We also perform the Kolmogorov-Smirnov test (K-S), seeKolmogoroff (1941); Smirnov (1948); 
Massey (1951), which does not require binning the data. The K-S test, as implemented by 
the FORTRAN subroutine KSONE in Press et al. (1992), finds the maximum distance, D, 
between the theoretical and the astronomical DF as well the significance level Pks , see 
formulas 14.3.5 and 14.3.9 in Press et al. (1992). Values of Pks > 0.1 assures that the fit is 
acceptable. 



5. Astrophysical applications 

This section reviews the galactic IMF as modeled by three and four power laws PDFs 
and fits the masses of the cluster NGC 2362 and the cluster NGG 6611 with the various 
PDFs here considered. 



5.1. Galactic IMF 

The IMF is usually modeled by two or three power laws of the type 

Pstars{m) OC X""* , (52) 

each zone being characterized by a different exponent CKj. In order to have a PDF normalized 
to unity, one must have 

/ Cim"^'dm^l . (53) 



i=l,3 
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For example, we start with ci—1: will be determined by the following equation 

ci(0.5-e)-"i =C2(0.5 + e)-"2 , (54) 

where e is a small number, e.g., e = 10"'^. In the previous equation wc insert a\ = 1.3 
and ^2 = 2.3 and therefore C2 = 0.503. The same procedure applied to C3 gives c^= 0.506. 
The integral of Pstarsi^) over the field of existence now gives 4.14, but according to the 
requirement of normalization as given by Equation (53), it should be 1. As a consequence, 
the three constants are now ci = 0.24, C2 = 0.1205, and C3 = 0.1206, which is the same as 
equation (59) in Kroupa et al. (2012) 

0.24X-1-3 if O.O7M0 < m < O.5M0 
p{m) = { 0.12X-2-3 if O.5M0 <m< I.OMq (55) 
0.12 x-^-^ if I.OM0 < m < IOM0 . 

The mean of the galactic IMF is given by a numerical integration over the three zones 

rh=^ Cimm~"*dm = O.389M0 . (56) 

The presence of the brown dwarfs means the use of four power laws instead of three power 
laws: 

'2.194 if O.OIM0 < m < O.O7M0 

, , 0.153 if O.O7M0 < m < O.5M0 , , 

p(m) = { U - U 57 

^ 0.076x-2=^ if O.5M0 < m < I.OMq 

_0.076x-2-^ if I.OM0 < m < IOM0 , 

where in order to have a continuous PDF, the BDs have the range O.OIM0 < m < 0.07Mq 
rather than O.OIM0 < m < O.15M0, see equation (59) in Kroupa et al. (2012). We have 
covered the galactic four power laws, we now introduce the generalized four power laws 
PoiiTi; —CKi, —a2, — eta, —a^, mi, m2, ma, m^, m^) which in the case of NGC 2362 is 

PG(m; -0.01, -0.02, -1.1, -2.7, 0.01, 0.07, 0.50, 1.0, 10.) A^GC 2362 case , (58) 

and in the case of NGC 6611 is 

PG(m; -0.01, -0.6, -2.4, -2.7, 0.01, 0.07, 0.50, 1.0, 10.) NGC mil case . (59) 



5.2. IMF of NGC 2362 

A photometric survey of NGC 2362 allows of deducing the mass of 271 stars in the range 
1.47M0 > M > O.11M0, see Irwin et al. (2008) and the data in J/MNRAS/384/675 at 
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the Centre de Donns astronomiques de Strasbourg (CDS). Table 1 shows the values of Xred^ 
the AIC, the probability Q, of the astrophysical fits and the results of the K-S test. 

Table 1: Numerical values of Xred^ AIC, probability Q, D, the maximum distance between 
theoretical and observed DF, and Pks , significance level, in the K-S test for the mass 
distribution of the NGC 2362 cluster data (272 stars). The number of linear bins, n, is 20. 



PDF 


parameters 


AIC 


^red 


Q 


D 


Pks 


lognormal 


(7=0.5, uln — —0.55 


37.64 


1.86 


0.013 


0.07305 


0.10486 


double 
Pareto-lognormal 


(7=0.44 ,hln = -0.52 
a=5, [3 = 5 


40.42 


2.02 


0.008 


0.066103 


0.17882 


general beta 


a = 0.12,6 =1.47 
a = 1.67,/3=2.77 


29.09 


1.31 


0.17 


0.059141 


0.288813 


general beta 
+normal (NB) 


a = 0.12,&=1.47 
a = 1.67,^=2.77, a = 0.001 


31.09 


1.40 


0.13 


0.06412 


0.20612 


left 

truncated beta 


= 0.12,6=1.47 
a = 2.23,^=3.09 


31.19 


1.44 


0.1 


0.06158 


0.24572 


four 
power laws 


Eqn. (58) 


77.608 


4.89 


1.1710-^ 


0.16941 


2.60363 10-^ 



Figure 1 shows the fit with the left truncated beta distribution of NGC 2362 and Figure 
2 visually compares the four types of fits for NGC 2362. 

5.3. IMF of NGC 6611 

The massive young cluster NGC 6611 has been carefully analyzed from the point view 
of the IMF in the range 1.5M0 > M > O.O2M0. This means that also the BD range 
is covered, see more details in Oliveira et al. (2009) with data in J/MNRAS/392/1034 at 
the CDS. Figure 3 shows the fit with the left truncated beta distribution of NGC 6611 and 
Figure 4 shows a visual comparison of four types of fits for NGC 6611. Table 2 shows the 
values of Xred^ AIC, and the probabihty Q of the astrophysical fits and the results of 
the K-S test. Figure 5 shows the residuals and ^ function of the middle value of the 

logarithmic bin considered, both for the left truncated beta and for the lognormal. 
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Table 2: Numerical values of X^ed' AIC, probability Q, D, the maximum distance between 
theoretical and observed DF, and Pks , significance level, in the K-S test for the mass 
distribution of NGC 6611 cluster data (207 stars + BDs). The number of linear bins, n, is 
20. 



PDF 


parameters 


AIC 


Xred 


Q 


D 


Pks 


lognormal 


(7=1.029,/iLjv = -1.258 


71.24 


3.73 


1.310-^ 


0.09366 


0.04959 


double 


a= 0.979 ,hln = -1.208 


70.3 


3.89 


2.1310-^ 


0.07995 


0.13523 


Pareto-logiioriiial 


(\ = [, 3 = 1 












general beta 


a = 0.019,6 =1.46 
a = 0.56,/3=1.55 


39.29 


1.956 


0.0123 


0.11456 


0.007924 


general beta 


a = 0.019.6 =1.46 


41.3 


2.08 


0.008 


0.09476 


0.04545 


+normal (NB) 


a = 0.56,/3=l. 55,(7 = 0.001 












left 


a = 0.019,6 =1.46 


42.09 


2.13 


0.005 


0.06839 


0.27781 


truncated beta 


a = 0.55,/3=1.6 












four 


Eqn. (59) 


81.39 


5.18 


2.4119-9 


0.12514 


2.7239 10-3 


power laws 















6. Conclusions 

Motivations In the last 50 years, the IMF has been modeled progressively by one 
power law, by two power laws, three power laws, and four power laws. The three power 
law distribution has seven parameters and the four power law has nine, and they both have 
a finite range of existence. A second widely used fitting function is the lognormal, which 
is characterized by two parameters and is defined on the interval [0,oo]. In this paper, we 
have described a left truncated beta PDF which has (i) a lower and an upper bound, (ii) 
a finite value of probability on the lower bound rather than zero, (iii) two parameters, a 
and P, which fix the shape of the distribution, (iv) an analytical expression for the average 
value. Two physical meanings are distinguished: (i) the upper limit of the left truncated 
beta is connected with the maximum stellar mass, which is ~ 6OM0, (ii) the lower limit 
is connected with an unknown physical mechanism which limits the distribution in masses. 
Further on we remember that the lognormal PDF has the important disadvantage of missing 
the well-accepted Salpeter-type high-mass power law. The high masses behavior of the 
various PDFs here analyzed is reported in Figure 6 for NGC 6611 where the Pareto and 
truncated Pareto PDFs are evaluated for M > 0.43Mq, which means Salpeter slope -2.3 . 
From the previous Figure the discrepancy of the lognormal and double Pareto-lognormal at 
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high masses is evident. 

Goodness of fit tests The statistical tests here performed are spht in two: (i) a first 
test requires binning the data in order to evaluate x^, and the indicators are Xred^ AIC, 
and the probability Q; (ii) the K-S test does not require binning the data, and the two 
indicators are D and Pks ■ These two tests, when applied to NGC2362 and NGC6611, 

indicate that the beta family (general and left truncated) performs better than the lognor- 
mal distribution both when the binning of the data is computed, see Tables example, the 
K-S test for the mass distribution of NGC6611 indicates a confidence level of 27% for the 
left truncated beta and 5% for the lognormal. New confidence levels can be found with 
the Anderson-Darling test which is a modification of the K-S test, see Stephens (1974) 
and the discussion at https://asaip.psu.edu/Articles/beware-the-kolmogorov-smirnov-test. 
Currently tables of critical values for the Anderson-Darling test are available for the log- 
normal PDF but the critical values for other PDFs here explored are not yet available, see 
http:/ /www.itl.nist.gov/div898/handbook . 

Convolution The random sum (convolution) of a general beta and a normal random 
variable, as represented by Equation (45), when apphed to NGC6611 introduces a further 
parameter, a, which increases xled from that of the general beta, but decreases 

D and Pks in the K-S test, see Table 2. 

The mode A careful attention should be paid to the falloff of the IMF towards the 
brown dwarfs. The left truncated beta PDF, see PDF (34), once the numbers of the open 
cluster NGC6611 are inserted, see Figure 3, decreases after the maximum at m f=:i 0.019Mq. 
This fact can be explained by the following Taylor expansion 



The previous decreasing function converts itself into an increasing function when the inte- 
gration is performed 



/t(x; 0.019, 1.46,0.715,2.185) 
2.65 - 21.64 {x - 0.039) + 336.35 {x - 0.039)^ + O {{x - 0.039)^) 



(60) 




0.019 



/t(x; 0.019, 1.46, 0.715, 2.185) 



3.50 x - 10.82 x'^ + 112.11 {x - 0.039)^ 



(61) 



and we recall that the evaluation of the frequencies corresponds to an integration. 



Lognormal family The recently formulated double Pareto-lognormal distribution draws 
attention to a possible alternative to the lognormal. Our tests show that the double Pareto- 
lognormal lowers the value of the maximum distance, D, of the K-S test, see Tables 1 and 



- 15 - 



2. Inconveniently, at the moment of writing there are no analytical evaluations of the four 
parameters which characterize the double Pareto-lognormal. 

The astronomical sample The new PDFS here presented can be tested on an astro- 
nomical sample representative of the IMF. Currently not all the various catalogs available 
on CDS report the column of the mass. As an example the promising IMF of IC 348 , see 
Figure 11 in Alves de Oliveira et al. (2012) , is not available on CDS. 
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Fig. 1. — Logarithmic histogram of mass distribution as given by NGC 2362 cluster data 
(272 stars) with a superposition of the left truncated beta distribution when the number of 
bins, n, is 12, a= 0.12, b= 1.47, a — 2.23 and (3 =3.09. Vertical and horizontal axes have 
logarithmic scales. 
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Fig. 2. — Histogram (step-diagram) of mass distribution as given by NGC 2362 cluster 
data (272 stars) with a superposition of the left truncated beta distribution (full line), the 
lognornal (dashed), the double Pareto lognormal (dotted) and the four power laws (dot- 
dash-dot-dash) . Vertical and horizontal axes have logarithmic scales. 
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Fig. 3. — Logarithmic histogram of mass distribution as given by NGC 6611 cluster data 
(207 stars + BDs) with a superposition of the left truncated beta distribution when the 
number of bins, n, is 12, a= 0.019, b= 1.46, a — 0.55 and /3 =1.6. Vertical and horizontal 
axes have logarithmic scales. 
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Fig. 4. — Histogram (step-diagram) of mass distribution as given by NGC 6611 cluster 
data (207 stars + BDs) with a superposition of the left truncated beta distribution (full 
line), the lognormal (dashed), the double Pareto lognormal (dotted) and the four power laws 
(dot-dash-dot-dash). Vertical and horizontal axes have logarithmic scales. 
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Fig. 5. — The residuals of the fits to NGC 6611 cluster data when 12 logarithmic bins are 
considered. The empty stars represent the left truncated beta PDF and the filled triangles 
the lognormal PDF. 
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Fig. 6. — Survival function of NGC 6611 cluster data as logio ~ ^oqiq plot when M > 
O.43M0 : data (empty circles), survival function of the truncated Pareto pdf (full line) 
(a=0.43,b=1.46,c=1.3) and survival function of the Pareto pdf (dashed line) (c=1.3 , Salpeter 
slope -2.3). The left truncated beta distribution (dot-dash-dot-dash) , the lognormal (dot- 
ted) and the Double Pareto-lognormal ( dash-dot-dot-dot) cover all the range in mass with 
parameters as in Table 2. 



